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Abstract 

We study the neutralino sector of the left-right supersymmetric model. In addition to the possi- 
bilities available in the minimal supersymmetric model, the neutralino states can be superpartners 
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\q \ of the U(1)b-l gauge boson, the neutral SU(2)r neutral gauge boson, or of the Higgs triplets. We 



■^j- ' analyze neutralino masses and determine the parameter regions for which the lightest neutralino 

o" 

^sO ■ can be one of the new pure states. We then calculate the relic density of the dark matter for each 
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(-H , of these states and impose the constraints coming from the p parameter, the anomalous magnetic 



moment of the muon, b — ► 57, as well as general supersymmetric mass bounds. The lightest neu- 
tralino can be the bino, or the right-wino, or the neutral triplet higgsino, all of which have different 
couplings to the standard model particles from the usual neutralinos. A light bino satisfies all the 
experimental constraints and would be the preferred dark matter candidate for light supersym- 
metric scalar masses, while the right-wino would be favored by intermediate supersymmetric mass 
scales. The neutral triplet Higgs fermion satisfies the experimental bounds only in a small region 
of the parameter space, for intermediate to heavy supersymmetric scalar masses. 
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I. INTRODUCTION 

Recent measurements from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite 
[l| have shown that the cold dark matter (CDM) abundance in the universe is VLcDkth 2 = 
0.1 126^009045 wnere &>cdm is the density of the CDM species versus the critical density, 
and h is the present value of the Hubble parameter. This corresponds to about 22% of the 
energy density of the universe being present in the form of dark matter. 

The question of its nature and composition has been open problem for some time. A 
candidate for dark matter must be stable, or long-living; it must be electrically and color 
neutral, as required by astrophysical constraints; and it must be non-relativistic, which 
means that it should be massive. Since the relic density can be connected to the (thermally 
averaged) annihilation cross section of a dark matter candidate, its value indicates that these 
particles are interacting weakly. That is, dark matter candidates must be weakly interacting 
massive particles. Scenarios beyond the standard model provide several such exotic particles. 
Supersymmetry, in particular, provides the lightest neutralino as the lightest supersymmetric 
particle (LSP), which must be stable if R-parity is conserved. 1 

The MSSM Lagrangian, based on the symmetry group SU(3)c <8> SU{2) L ® U(l)y is 
invariant under the R-parity discrete symmetry, under which all standard model particles 
are even, and their superpartners are odd. This parity is an essential element of the MSSM, 
and forbids baryon and lepton number violating renormalizable couplings in the superpo- 
tential, which would lead to fast proton decay. Thus the R-parity is natural, and so is the 
possibility of having an LSP which emerges as a natural candidate for dark matter. In most 
supersymmetric scenarios, the LSP is the lightest neutralino 

Unfortunately there are problems with this CDM assignment in MSSM |3|. If one calcu- 
lates the relic abundance of the lightest neutralino assumed to be a pure bino (the super- 
partner of the gauge boson of the U(l)y gauge group), the relic density is too large (see 4( 
and references therein for a detailed discussion). To avoid this problem, one considers that 
there are several other supersymmetric particles of masses close to the LSP. If such sparticles 
exist, their co-annihilation with the LSP leads to a reduction of the LSP relic abundance 
Q. Extent studies of such phe_ exist, and several publicly available codes such as 



1 Gravitinos and Axinos have also been studies as dark matter candidates, but we shall not consider this 
possibility here. 



DarkSUSY p and mi cr OMEGA |7| , which include all the relevant coannihilation channels, are 
available for calculations within the MSSM. Cosmological and phenomenological constraints 
(such as b — > 37 J8|, (g — 2) M J9j], and Higgs mass bounds |13|) are also included. We could 
of course also abandon the pure bino state, and thus the constrained MSSM (CMSSM) 
which predicts it, and explore other versions of the lightest neutralino. However, this does 
not solve the problem because if one allows the lightest neutralino to be mostly left-wino, 
or higgsino, the relic density becomes too small. Ways to avoid this exist as well, such as 
considering non-universal gaugino masses, abandonin g th e pure states and examining mixed 
states (such as bino-higgsino) , or more annihilations |lfl ]. 

An alternative option is to explore scenarios beyond the MSSM. The extensions of the 
MSSM, motivated by the need to a dynamical solution to the naturalness problem associated 
with the /i parameter, offer novel CDM candidates not necessarily be long ing to gaugino or 
Higgsino sectors of the MSSM. In next-to-minimal model (NMSSM) 1J|, for instance, the 
neutralino can be extremely light because it is dominated by its singlino component, which 
does not couple to SM particles (except for Higgs doublets). The existence of a very light 
CP-odd Higgs boson provides the possibility for this neutralino to annihilate sufficiently to 
avoid b eing overproduced in the early universe. Similar features are also found in U(l)' 
models 11 1 or in more general extensions [12j. This raises the hope that models beyond 
MSSM, preferably comprising the neutrino masses as well, can provide viable scenarios for 
alternative neutralinos which have different features than those in MSSM. 

We propose here to look at the left-right supersymmetric model (LRSUSY) for candidates 
for dark matter. The advantages of the LRSUSY are many, such as combining supersymme- 
try with right handed symmetry, and thus providing a mechanism for neutrino masses with 
the seesaw within a supersymmetric scenario. From the point of view of dark matter, the 
LRSUSY model is interesting because, unlike in MSSM, the usual explicit R-parity violating 



terms are forbidden in the superpotential by the symmetries of the model 
model naturally predicts a lightest supersymmetrical particle. 
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2 Thus the 



The LRSUSY model emerges naturally in some superstring theories |17| . or in the break- 
ing of such SUSY GUT scenarios, such as SO(10) |l8[. LRSUSY is based on the symmetry 
group SU(3) (g) SU{2) L ® SU{2) R <S> U(1)v=b~l an d thus provides several new sources for 



2 Note however that the R-parity can be broken spontaneously [16} ■ 



exotic neutralinos, in the (neutral) partners of the SU(2) R or U(1)b-l gauge bosons, as 
well as the fermionic partners of several Higgs bosons. We analyze the candidates for dark 
matter in this model to see if they can avoid the problems that plague MSSM. After briefly 
reviewing the model in Section [HJ we perform a comprehensive study of neutralino mass 
eigenstates in Section 11111 followed by a calculation of relic density in Section I1VI We look 
at the simplest scenario, that of pure states and avoid co-annihilation with other supersym- 
metric particles by choosing the supersymmetric masses accordingly. We include constraints 
from the WMAP, b — > 57, Ap, and muon g — 2 as well as from experimental bounds on 
supersymmetric masses from direct collider searches. 

II. THE LRSUSY MODEL 

The most general superpotential for the group SU(3) <S> SU(2)l <S> SU{2) r ® U{1)v=b-l is 
Q: 

W = Y®Q T %iT 2 Q c + YfL T ^ir 2 L c + i(Y lr L t t 2 5 l L + Y LR L cT r 2 A R L c ) 

+fi LR [Ti(A L 5 L + A R S R )] + m/rrfaQTiTzQj) + W NR , (2.1) 

where Wnr denotes (possible) non-renormalizable terms arising from higher scale physics 



or Planck scale effects 
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2JJ. Here the matter fields are defined as 



Q = I j -(3,2,1,1/3), Q c = I d c I ~ (3*, 1,2, -1/3), 

~ (1,2, 1,-1), L c =( 6 \ -(1,1,2,1), (2.2) 

where the numbers in the brackets denote the quantum numbers under SU(3)c <8> SU{2)l ® 
SU{2) r ®U{1)b-l- The Higgs sector consists of the bidoublet and triplet Higgs superfields: 



1,2,2,0), $ 5 



1,3,1,-2), ^==| f L 

(1,1,3,-2), 5 R =~- [ f R . + 1 -(1,1,3,2). (2.3) 

°R ~V2 







where Al and Sr are introduced in the model to cancel the fermionic anomalies introduced 
by the fermionic partners of 5l and Ar. The vev's of the Higgs fields in the LRSUSY can 
be chosen 
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h)= }, (A fl > = Ar }, (S R ) = . (2.4) 

Some comments and explanations about the vev's chosen are required: K\ and k 2 are the 
vev's of the MSSM-like Higgs bosons. They are responsible for giving masses to the quarks 
and leptons and they also contribute to M\y L - We take the vev's of cj)\ 2 and 0°i to be zero 
because they induce FCNC at tree level (in both the leptonic and hadronic systems), as 
well as being responsible for Wl — Wr mixing. The vev's could also have a phase which 
induces CP violation, which is severely restricted in the kaon system. The non-MSSM Higgs 




vev's, yx T and va r are responsible for neutrino masses. For one generation (see J22J, also 

HflQ) 

(Y {u) ) 2 k 2 
m v = Y lr v Sl - K L ' 1 (2.5) 

where v$ L must be very small and va r large, phenomenologically . In addition, va l and 
v$ L enter in the formula for the mass of Wl (or the p parameter), while Va r ,Vs r enter in 
the formula for the mass of Wr. It is thus justified to take v/\ r , vj ,j to be negligibly small. 
For va r there are two possibilities: either va r is ~ 10 13 GeV J2J, |25J], which supports the 
seesaw mechanism, leptogenesis and provides masses for the light neutrinos in agreement 
with experimental constraints, but offers no hope to see right-handed particles; or va r is 
pa 1 — 10 TeV, but one must introduce somethin g el se ( g enerally an intermediate scale, or 
an extra symmetry) to make the neutrinos light |l5J . I20l l25| . Note that both va r and v$ R 
contribute to the mass of Wr 2j| , but only one needs to be heavy. Since va r is responsible 
for heavy right-handed neutrino masses, it must be large. Thus v$ R is not constrained by 
the data. The mass term for neutralinos is given by 



C N = 



l -i)° T Yi)° + H.c. , (2.6) 



where ^° = {-i\ v , -i\° L , 4>° u , 0° 2 , A% 5° R , -i\% 0? 2 , 4>° 2l ) T . Here \ v is the U(1) B _ L bino, \° L 
the left-handed neutral (SU(2) L ) wino, and A^ the right-handed neutral (SU(2) R ) wino. 



The rest of the fields are higgsinos. The neutralino mass matrix Y is equal to j2 
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where /% = /ii, (i 7^ j) and \ilr = ^2 are assumed. 3 The upper 4x4 part of the matrix 
contains MSSM-like states. There are still too many unknown parameters and for further 
simplification one can define the followings: 



9l = 9r = 9-, 9v 



KT 



k\ = 174 GeV, tan/5 = — . (2.8) 



Vcos 26 w 

If we assume that GUT relations between gaugino masses hold, we can simplify the param- 
eters further, but for now we keep them general and discuss specific scenarios later. 

Two roots of the characteristic equation for Y are already known, ±/ii, and we are left 
with a 7 th order equation to solve. There is no exact solution and a numerical or semi-exact 
approach is necessary. Whereas there are constraints on Mw R and Mz R |28|, there are no 
constraints on My, the U(1)b-l gaugino mass, or on M L ^, the SU(2) L ^ gaugino mass 
parameter. 

Basically the number of free parameters can be chosen as My, Ml, M r , va, vs, /ii, fi2, 
and tan/3. So, we could look at under what circumstances the lightest neutralino becomes 
mainly pure state (bino, left-wino, right-wino, or higgsino) or a mixed state as bino-higgsino 
and right- wino-higgsino. 



From now on we drop the subscript "R" from va r , vg R . 



III. THE CLASSIFICATION OF SCENARIOS 

Starting from the neutralino mass matrix, we can classify the following cases: 4 

1. The lightest neutralino is mostly bino (Xy in our notation). Note that this the the 
B — L gaugino, not to be confused with the U(l)y MSSM bino. To obtain this bino 
as the lightest state, we must have My « 50—100 GeV and v$ ~ 50 — 100 GeV 
(both light) while \ii ~ 1 TeV or larger (heavy). Decreasing /12 to 200 — 500 GeV, the 
lightest state becomes a mixture of Xy and 8^. 

2. The lightest neutralino is mostly right-handed wino (X° R in our notation). To obtain 
this, the masses of the other two gauginos must be larger than that of X R : My, Ml rs 
600 TeV and v s « 50-100 GeV is light; while \i x « 1 TeV, /x 2 » 3-5 TeV. Decreasing 
/i 2 to 1 TeV, the lightest state becomes a mixture of X R , Xy and 8 R . 

3. The lightest neutralino is the right-handed higgsino 8 R . This scenario is obtained for 
a large range of parameters, as long as both \xi and v$ are small, ~ 200 GeV or larger, 
and (0, 100) GeV, respectively. In fact, this requirement is satisfied for a wide portion 
of the parameter space, as long as \x^ and vg are smaller than other parameters in the 
theory. This scenario is interesting since 8 R does not couple to any Standard Model 
particles (or their SUSY partners). As explained previously, Vg is not constrained by 
the data. 

4. The lightest neutralino could be a mixture of Xy, 8 R and A° R . In this case we can 
get eigenvectors and eigenvalues analytically and calculate the relic density for the 
mixed state. This is a generalization of the first case, where both A R and 8 R are 
included. Here we can take advantageous ratios of the vevs of 8^ and A° R to get only 
one combination of 8% and A° R mixed with the bino. This is the case if we assume the 
vev's vs and v& equal. 

5. The lightest neutralino could be a mixture of A^ (the right-handed wino), and the 
non-MSSM higgsinos 8^, A° R . As opposed to the B — L bino, the right-handed wino 
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We set i'a = 1.5 TeV for all the cases unless otherwise stated. Note also that there will be no rigorous 
model constraint analysis on the masses of the LSP in this section. We leave that to the next Section. 



mixes with MSSM higgsinos as well, which will be considered in a separate scenario. 
Again like in scenario (4), one can make the right- wino couple with one combination 
of S R and A R if their vev's are assumed to be equal. 

6. Finally, the lightest neutralino could be a mixture of X R with only the MSSM higgsinos, 
(/>5i and </> 22 . This can be the case if both Mr and [i\ are small compared to Ml, My, fi2, 
and the vev's Va and v§- One could consider v^ = v$ case and decouple 5 R and A R 
by taking /x 2 to be large. 

We are not interested in the MSSM lightest neutralino scenarios (in which the left-wino 
or MSSM higgsino mixed with left wino are the LSP), since these have already been studied. 
As far as we can tell, scenarios 1-6 are the important (most striking) possibilities. 

A. The lightest neutralino state 

The following four sets of figures further illustrate the scenarios mentioned above. We 
assume the composition of the lightest state as (for the first three Scenarios) 

\Xi) = Nu\X°l) + N 12 \X° R ) + JV 13 |Av) + N u \4>° u ) + N 15 \4>° 22 ) + A 16 | A ^) + N l7 \~5° R ). (3.1) 

First we look at the possibilities for Scenarios (1) and (3). Here we want the bino to be 
the lightest, so we take Ml, Mr and Hi large. Accordingly, we need Vg to be small. Varying 
fj,2 will take us from a mostly bino lightest state to a mostly S R lightest state. 

In the left panel of Fig. Q the difference between the bino and higgsino compositions 
of the lightest neutralino state is shown as a function of /i 2 by choosing (My,v$) = 
(0,0), (0, 100), (0,400), and (200,400) GeV. The rest of the parameters are fixed as M L = 
Mr = 600 GeV, /ii = 500 GeV, v& = 1.5 TeV, and tan/? = 2. In the right panel, the mass 
of the lightest neutralino is given as a function of /i 2 for the same parameter values. As seen 
from the figure the lightest state is pure higgsino 5 R for very small /i 2 values, regardless of 
the values of (My, v$), as long as they are smaller than 500 GeV. The ratio is more sensitive 
to the vev of the right-higgsino (v$) than to the U(1)b-l gaugino paramater My. At large 
/i 2 the state becomes pure bino from pure higgsino and its mass shows a strong dependence 
on /i 2 . So, one can consider these two limiting cases as realizations of the scenarios (1) and 
(3), mentioned above. From the right panel, one can conclude that, except for very small 

8 



values of /i2, the cases where My and v$ are larger than 200 GeV predict a lightest neutralino 
with a mass in the range of 200-300 GeV. Otherwise its mass remains less than 150 GeV. 
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FIG. 1: On the left panel, the difference between bino and higgsino compositions of the lightest 
neutralino state as a function of [12 for various (My,v$) values in TeV, for Ml = Mg = 600 
GeV, /ii = 500 GeV, v& = 1.5 TeV, and tan/3 = 2. On the right panel, the mass of the lightest 
neutralino as a function of fi2 for the same parameter values. 



Next we look at Scenarios (2) and (3). Here we want the right-wino to be the lightest 
neutralino, so we take Ml, My and [i\ large, and Vg to be small. Varying ^2 will take us 
from a mostly right-wino lightest state to a mostly 5 R lightest state. 

In Fig. |2J we plot the difference between the bino and the higgsino compositions of the 
lightest neutralino state as a function of \ii for various (Mr, v$) values, for Ml = My = 600 
GeV, Hi = 500 GeV, v& = 1.5 TeV, and tan/3 = 2. On the right panel, the mass of the 
lightest neutralino is shown as a function of fi 2 for the same parameter values. The third 
diagram shows the bino composition of the lightest state. This case exhibits a very similar 
dependence on /i 2 as the bino-higgsino case. The curves for various (Mr, v$) pair indicate 
that the ratio is insensitive to their chosen values, as long as they are assumed to be less 
than 200 GeV. The bino composition becomes significant only around /i2 = 1 TeV, where 
the right-wino and higgsino mix almost equally, and it is negligible as /12 becomes larger. 
The contribution of the bino with respect to that of the right-wino is significant for very 
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FIG. 2: On the left panel, the difference between right-wino and higgsino compositions of the lightest 
neutralino state as a function of fi2 for various (Mr,vs) values in TeV, for Ml = My = 600 GeV, 
fi\ = 500 GeV, va = 1-5 TeV, and tan (5 = 2. On the right panel, the mass of the lightest neutralino 
as a function of /i2 for the same parameter values. The third diagram shows the difference between 
bino and higgsino compositions of the the lightest state. 

small fi 2 values (~ GeV) but such small values are excluded from the lower bound on 
the mass of the lightest chargino, which is around 90 TeV. The mass of the lightest state is 
less than 150 GeV for v 5 < 200 GeV if M R is very small, or for (M R , v s ) < (50, 50) GeV. 
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Otherwise, it is larger than 200 GeV for \x-i larger than 500 GeV. While the composition 
of the state is insensitive to Mr and Vg, the mass of the lightest neutralino is sensitive to 
both. Here one can recover Scenarios (2) and (3) in the limiting cases and the mixed state 
(\X° R ) + \5° R )) in between. 

Finally, we look at realizations of Scenarios (4), (5), and (6). Here we are interested in 
mixtures of bino-(non MSSM) higgsino, right-wino- (non MSSM) higgsino, and right-wino- 
MSSM higgsino. We need to decouple the MSSM particles for scenarios (4) and (5), so 
we take Mi and [X\ to be large. In addition we take V& = v§ (as opposed to the previous 
scenarios, where we took one of them small). We vary My, Mr and [i-i to go from one case 
to another. 

The realization of the last three scenarios discussed above requires some radical changes 
in the parameter set and has a distinct structure. We first decouple the ^>\ 2 and <\P 2 \ fields. 
Then, without loss of generality, we rotate the basis 

" A°>, |A£), |A v ), |ft>, |&>, lA^), \5 R }\ - (|A°), \X r ), \\ v ), |#.), \$L), \ffi), m} (3.2) 



such that 



:) = 



I0?i) ± |&> 



fT±) = 



\&h) ± 14.) 



V2 ' " ±/_ V2 

In this new basis the mass matrix Y in Eq. (|2.7|) becomes 
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(3.4) 



where we have assumed va — v$ = vr, and further defined k = K\ ± k^- Under this 
assumption the rotated state ff, decouples and only fj°_ remains, which allows us to analyze 
scenarios (4), (5) and (6) more conveniently. 

Fig. E3 shows the difference between the right-wino (bino) and the rotated non-MSSM 
higgsino (fj°_) compositions of the lightest neutralino state as functions of /i2 by choosing 
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FIG. 3: On the left 'panel, the difference between right-wino and rotated non-MSSM higgsino com- 
positions of the lightest neutralino state as a function of \ii for various {Mr, My) values in TeV, 
for Ml = 600 GeV, ix\ = 5 TeV, va = vs = 1 TeV, and tan/3 = 2. On the right panel, the 
difference between the bino and the rotated non-MSSM higgsino compositions as functions of [12 
for the same parameter values. The third diagram shows the mass of the lightest neutralino as a 
function of [i<i . 

various (Mr, My) values in the left (right) panel. The rest of the parameters are taken 
as M L = 600 GeV, (n = 5 TeV, v A = v s = v R = 1 TeV, and tan/5 = 2. Here l^] 2 
is the amount of the rotated higgsino field fj°_ in the lightest neutralino state. One can 
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divide the discussion into two regions-one with small Mr < 50GeV and My < 1 GeV, and 
the other with larger Mr and My. In the first part, neither the right- wino nor the bino 
is the lightest neutralino till /x 2 ~ 500 GeV. The decoupled state 77° is the lightest state 
(which can not be seen from the figures 5 ). For fi 2 > 500 GeV, the lightest state is mainly a 
right- wino or a right- wino-higgsino mixture. Otherwise, for intermediate values of fi 2 , it is a 
mixture of right-wino, bino, and higgsino 77° . For the second part where Mr > 50 GeV and 
My > 1000 GeV, the curves are horizontal lines passing through zero for both graphs until 
\i 2 ~ 2 TeV. In that region of the parameter space the lightest state is pure left-wino. The 
right- wino-higgsino mixture is possible as the lightest neutralino state only after that point. 
The bino-higgsino lightest state is not realized in this region. From the third diagram in 
Fig. |3J one can say that within the range considered for /j, 2 , the mass of the lightest state is 
greater than 300 GeV for /i 2 > 500 GeV. So, the scenarios (4) and (5) require the lightest 
state to be rather heavy unless fi 2 is too large. This takes us to pure right-wino scenario. 
Including all the possible model constraints in the next section will restrict our parameter 
space further. 

In Fig.HJ the difference between the composition of the lightest neutralino as a right-wino, 
or a MSSM-higgsino is shown on the left panel as a function of Mr for various fi 1 values in 
the range to 500 GeV. The other parameters are chosen as Ml = My = 1 TeV, fi 2 = 20 
TeV, Va — v§ — 1 TeV, and tan/3 = 2. The ratio is sensitive to both Mr and /ii. For values 
of \i\ smaller than 25 GeV, the lightest state is pure MSSM-higgsino and for values larger 
than 450 GeV the state almost is pure right-wino, if Mr is equal to, or less than, 150 GeV. 
Evenly mixed states are obtained for intermediate values of \i\. As Mr becomes equal to, 
or larger than, 400 GeV, the lightest state remains pure right-wino for the most part of the 
/ii range. As suggested by the right panel, the mass of the lightest state is very sensitive to 
the parameter Mr and can be quite large for large Mr and \i\ values. 

This completes analyses of possible LSP candidates in LRSUSY for certain patterns of 
the model parameters. The figures illustrate nature and purity of the lightest neutralino 
states as a function of various gaugino masses and Higgs vevs. In the next section we will 
perform a detailed analysis of the relic abundance of LSP and its ability to explain the CDM 



5 Note that the horizontal line passing through zero does not always mean equal mixing. One could obtain 
such a result for vanishing individual contributions as well, which is indeed the case here. 
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FIG. 4: On the left panel, the difference between right-wino and rotated MSSM higgsino com- 
positions of the lightest neutralino state as functions of /ii at various Mr values in GeV, for 
Ml = My = 1 TeV, /i2 = 20 TeV, v/_ = v$ = 1 TeV, and tan/3 = 2. On the right panel, the mass 
of the lightest neutralino as a function of fi\ for the same parameter values. 

in the universe. 

IV. RELIC ABUNDANCE ANALYSIS OF CDM 

In this section we calculate the relic density of cold dark matter within the LRSUSY 
framework. Before presenting our results we give some details of the procedure followed. 

The time evolution of the number density rii for a relic particle is given by the Boltz- 
mann equation (taken within spatially flat Friedman- Robertson- Walker background). Fur- 
thermore, a single Boltzmann equation can be defined J5| for the total number density, 
n = J2i n i 



h + 3hn = - (cr e ffv) [n 2 - n^ q ] 



(4.1) 



where h is the Hubble parameter, v is the relative velocity of the annihilating particles, 
and n eq is the number density corresponding to the sum of each species number density at 
thermal equilibrium (that is, the density in the early universe). Here a e ff is the properly 
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averaged effective cross section of the CDM candidate into ordinary particles (i.e., SM par- 
ticles including the Higgs bosons) and will be defined below. Clearly, "( )" stands for the 
thermal average. 

In this study we concentrate on regions of the parameter space where the coannihilation 
effects are not si gnif icant. The relativistic thermally- averaged cross section times relative 
velocity reads as J29J 

/w dsa(s) y/s(s/± - m 2 ) K x (^/(m x x)) 
(a eff v)(x) = * — , (4.2) 

where m x is the mass of the annihilating particle, the lightest neutralino in our consideration; 
K\ and K 2 are modified Bessel functions; x = T jm x is the dimensionless temperature 
parameter; \fs is the center of mass energy of the annihilating neutralino pair; and a(s) is 
the cross section for the annihilation reaction x?Xi ~* ^sm, where Xsm represents all the 
allowed two-body SM final states. We will not approximate Eq. (I4.2J1 by expanding in v in a 
Taylor series, as this approximation might not always be safe, especially for cases with fast 
varying integrands. 

In general, there is no analytical solution to the Boltzmann equation, Eq. (|4.1|) . which is a 
Riccati-type equation. Thus a numerical approach is required. There exist several different 



.7 



ways to proceed in the literature |2l, l2i|. Here we define a freeze-out temperature parameter 
Xf = Tp/m x , which we use as an approximate solution to Eq. ()4.1|) |30j ] as 



m x Mpi\ 45 V T 7/ J X0 dx{cr e ffV)(x) 

where Mpi is the Planck mass, T is today's temperature (~ 2.742K°), so that x can be 
approximated as zero; T X (T 7 ) is the present temperature of the neutralino (Cosmic Mi- 
crowave Background), and g* represents the SM effective number of degrees of freedom at 
the freeze-out temperature (g* ~ 81 is used). Note that using an approximate solution of 
the Boltzmann equation instead of solving it numerically will introduce an uncertainity of 
up to 10% into our results. 

The freeze-out temperature parameter xp can be obtained from the following transcen- 
dental equation 

~m x Mpi j 45 



— = log 

xf 



2tt 3 
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v^^^^ 1 ^ 



(4.4) 



which can be solved iteratively. As a starting point we need to choose a value for xf- We 
used 1/25 in our calculations. The range for xf changes between 1/25 to 1/15. Finally, the 
relic density of neutralinos at present time is defined as 

l2 _ P( T o) 



tt x h< 



(4.5) 



where p(T ) = m x n x (T Q ) is the density of the neutralino, and p c = 3h 2 M Pl /87i = 8.098 x 
lCT 47 GeV 4 is the critical density of the universe. On the left hand side of the equation, h 
is the normalized Hubble expansion rate and its value today is 0.71. Thus from the central 



value of fi x h 2 , Q x is found to be about 22%. 
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FIG. 5: The annihilation Feynman diagrams contributing x?X? — * ^SM in the LRSUSY model. 
The MSSM contributions and u— channel diagrams are not shown here. Contributions are given in 
the basis of pure states whose mixtures can be easily constructed. Here f represent all fermions, 
and f all sfermions, h° is the SM-like Higgs particle, and 4>\,$_ are the rotated higgsinos fields 
defined in the previous section. 

To obtain our final results, a three dimensional integration needs to be carried out nu- 
merically. The three parts are: an integration over the Mendelstam variable t to compute 
a(s) for each subprocess in the kinematically allowed region, an integration over the center 
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of mass energy squared s from the threshold to practically infinity, and finally an integration 
over x to compute the thermal averaging from the freeze-out temperature to today's (xq), 
which we approximate to zero. Mathematica is used for the computation and the matrix 
element calculations have been carried out with FeynCalc 31 1. 

The Feynman diagrams contributing to the annihilation cross section are given in Fig. |3J 
We analyze here pure state contributions only, while any mixed case scenario can be cal- 
culated in a straightforward manner. Since our analysis concentrates on the non-MSSM 
scenarios, the MSSM contributions are not shown here. Depending of the center of mass 
energy available, the resonance problem in s-channel is handled using the Briet-Wigner 
prescription. 

The mass of the Z R boson in diagram (a) is M Zr = =(v$ + ^a) [26]. In 

V cos 28w 
diagram (b) we sum over all left- and right-handed quarks and leptons in the final state. 

The sfermion / has a (mass) 2 given by, for U-type squarks and sneutrinos: 

m s 2 usy + Ml {Tl - Q u sin 2 9 W ) cos 2(3 m Uk (A -ft cot (3) 



K4l DUD I ZJ L \ U -wu,- ,, , U k 

m Uh ( A - fj, cot p) M^sy + M 2 Zl Q u sin 2 9 W cos 2f3 



and for D-type squarks and sleptons 

m s 2 usy + M z L ( T d ~ Qd sin 2 9w) cos 2(3 m dk (A -ft tan (3) 



(4.6) 



r2 i\/r2 n.™2, 



M d Dk 

\ m dk {/i — /i tan p) iW susY + M ZiS d sm ~ ^ w cos ' Z P J 

(4.7) 



Here Msusy represents the universal scalar mass, A and \i the trilinear and bilinear scalar 
parameters, respectively, m u ^ are quark masses and the index k labels generations. It 
is customary to restrict the supersymmetric parameter space by drawing contours in the 
Mi/2 — ^susy plane, where M1/2 is the relevant gaugino or higgsino mass parameter and 
Msusy the relevant scalar mass. When the bino or right-wino are the LSP, we take their 
masses to be the significant gaugino mass. In diagram (c) the final state h° is the MSSM-like 
lightest Higgs boson, which we constrain to have a mass M^o = 115 GeV. In LRSUSY, its 
(mass) 2 is given by the lowest eigenvalue of the matrix J26J: 

where v 2 = vj + v\ + k\ + k 2 . and m^^ is the Higgs mass parameter. Note that while the 
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cross sections for the decay of Xy, X° R depend on Msusy, the one for the 5 R higgsino does 
not; we must find another relevant parameter for that case. 

V. NUMERICAL ANALYSIS 

Based on our classification of some possible mass scenarios in Section IHIi we analyze 
the first three mass scenarios with the lightest neutralino as a pure bino, right-wino, or 
higgsino state. The feasibility of the last three scenarios, assuming mixed gaugino-higgsino 
states, can be analyzed in the light of the pure state predictions. Throughout the numerical 
calculations, we have chosen the parameters of the model such that the Ap bound, taken as 
Ap < 0.002 [28], is always satisfied. 

Fig.lHlshows the relic density Q x h 2 of the neuralino as a pure bino state in the My — MgusY 

plane for fii = 500 GeV, A = 200 GeV, and tan/3 = 2. There is no s-channel contribution 

- f - 
and the t-channel fermion pair production XyXy — > ff, (/ = u,c,t,d,b,s...) is the only 

channel available. 6 Due to the s-wave supression, the fermion pair production for heavy 

fermions is dominant, but we include all possible channels. Unlike the MSSM bino, the 

B — L bino does not couple to W + W~ pairs or to Higgs pairs. For the exchanged sfermions, 

a flavor conserving scenario with a common scale Msusy is assumed. In the figure, the upper 

triangular region in green is excluded by the lower mass bound requirement on the lightest 

neutralino from LEP2 and by the lower bound of the lightest sfermion mass determined 

mainly by Msusy- (We calculated the masses of the sfermions and require that they are 

at least 15 GeV higher than the mass of the bino to avoid co-annihilation channels). The 

light blue regions are excluded by the WMAP measurements. In addition, the gray vertical 

strip on the left is excluded by the branching ratio of b — > 57 upper bound, for which the 

formulas are given in the Appendix. We assumed that 2.4 x 10 -4 < BR(b — > 57) < 4.1 x 10 -4 

range is allowed. In numerical calculation, we used the following set of parameter values, 

consistent with pure bino scenario: Ml = Mr = 3 TeV, rrig = 400 GeV, A = 200 GeV, 

pi = 500 GeV, /i2 = 4 TeV, v$ = 100 GeV, and v& = 1.5 TeV. 7 As seen from the figure 

BR{b — > sj) constraint excludes light M S usy (< 150 GeV) region. We also take into 



6 We recall that the bino here is the fermionic partner of the U(1)b-l gauge boson and is different from 
the bino in the MSSM. 

7 We ensure at all times that the bino is not just the lightest neutralino, but a pure state. 
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FIG. 6: The relic density of a pure bino state in My — Msusy plane for u\ = 500 GeV, A = 200 
GeV, and tan/3 = 2, in the scenario in which the bino is the lightest neutralino. The upper 
triangular region in green is excluded by the lower bound requirement of the lightest neutralino 
mass from LEP2 and by the lower bound on the lightest sfermion mass constraint. The light blue 
region is excluded by WMAP measurements. The gray vertical strip on the left is excluded by 
b — > sj. The white strip agrees with the WMAP at 2a level. The anomalous magnetic moment of 
the muon is satisfied within the entire region. 



account of the anomalous magnetic moment of muon a M = (g — 2)^/2, where the analytical 
formulas are given in the Appendix. The dominant chargino-sneutrino and neutralino- 
smuon loop contributions are included. In numerical computation, like BR(b — ► 57), the 
same parameter set is assumed and the range for the deviation from the SM, —0.53 x 10~ 10 < 
Aa M < 44.7 x 10 -10 , at 95% CL is used based on e + e~ data (see J35( and references therein). 
Aa M depends mainly on the chosen values of A, fj,i, tan/3, and Msusy- We observe that 
it remains in the allowed range throughout the entire (Msusy, My) interval considered, as 
long as Hi tan (3 — A > holds. 

Consequently, the white strip is the only region satisfying the WMAP at 2er level. A 
realization of such a scenario requires light sfermion masses in the range ~ 150 — 240 GeV 
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with an acceptable LSP mass in the region ~ 60 — 220 GeV. In the light blue region on the 
right, the relic density Q x h 2 \ x= \ v becomes larger than experimentally allowed values, while 
in the left painted region it is smaller than experimental requirements. This scenario can 
only be considered feasible in frameworks where at least there is at least one light sfermion, 
like stop t, or stau f. 

350 




200 300 400 500 600 700 800 900 1000 
M S USY [GeV] 

FIG. 7: The relic density of a pure right-wino state in Mr — MgusY plane for [i\ = 1 TeV, A = 200 
GeV, and tan/? = 2 in the scenario in which the right-wino is the lightest neutralino. The left upper 
small triangular and the lower regtangular regions in green are excluded by the lower bound mass 
requirement on the lightest neutralino from LEP2 and by the lower bound on the lightest sfermion 
mass. The light blue region is excluded by WMAP measurements. The grey region is excluded by 
b — > sj. The white strip agrees with the WMAP at 2a level. The anomalous magnetic moment of 
the muon is satisfied within the entire region. 



Next, in Fig.[7|we show the relic density £l x h 2 of the neutralino as a pure right-wino state 
in the Mr — M^usy plane for fii = 1 TeV, A = 200 GeV, and tan/3 = 2. The right-wino 
as a neutralino is a completely new possibility peculiar to the LRSUSY model. There is no 
s-channel contribution and the possible annihilation contributions are X R X° R —*■ ff, (/ = 



u, c,t, d, b, s. 



and X R X° R 



h h , when kinematically allowed. An important feature in 
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this case is that there exists an additional condition - the LEP lower bound for the lightest 
chargino, which we take as 90 GeV. We also insure that the lightest chargino remains 
always heavier than the lightest neutralino, at least 15 GeV heavier, to avoid significant 
contributions from co-annihilation channels. For the mixing in the chargino sector, we 
assume all parameters (except Mr and v$) large enough to obtain a pure right-wino state. 
Using vs = 30 GeV, a 100 GeV or larger mass is obtained for the right-wino, while fulfilling 
the bound on the mass of the lightest chargino. In order to obtain a neutralino lighter than 
the lightest chargino, M R should be at least 160 GeV or larger; hence the excluded region 
(in green) in Fig. For illustrative purposes only, we have included, in Table 1, the values 
of the masses of the LSP, as well as the masses of the second lightest neutralino, the lightest 
chargino, and the lightest sfermion, for all three scenarios investigated. As in the case of 
pure bino, the mass upper bound comes from sfermion sector. After implementing WMAP 
constraint, the only surviving region at 2<r level is shown in white strip. This scenario 
satisfies the WMAP conditions for a relatively large Msusy scale, 600 GeV or larger. Note 
that the upper bound for the mass of the right wino is closely related to the chosen values 
for the other gaugino masses, the vev's, and /12, assumed heavy. For larger Mr values, the 
lightest neutralino will not longer be a pure state. But one can see from Fig. [7| that it is not 
difficult to satisfy the WMAP result. Here too the constraints b — > 57 and (g — 2) of the 
muon are taken into account with the parameter set, Ml = My = 1.5 TeV, m~ g = 400 GeV, 
A = 200 GeV, \i x = 500 GeV, ll 2 = 4 TeV, v s = 100 GeV, and v A = 1.5 TeV. The b -> s 7 
constraint for this case excludes mainly Msusy < 600 GeV region but practically doesn't 
effect the WMAP allowed region, (g — 2)^ does not constrain the parameter space, as a light 
right-wino does not contribute significantly to the muon anomalous magnetic moment, and 
the other neutralinos are heavy. 

In the third scenario we investigate the possibility that the S R is the lightest neutralino. 
This an interesting scenario, as 5 R is only introduced to cancel anomalies in the fermionic 
sector, and, because of its B — L quantum number its direct interactions with matter are 
forbidden. The relic density prediction for such scenario is shown in Fig. |H1 as a function of 
the higgsino mass Mp , for v§ = 100 GeV, Va — 1-5 TeV. Only the s-channel contribution 

R 

is available for the annihilation process S R S R -5 //, (/ = u, c, t, d, b, s...). Note that for this 
process, the cross section and thus the relic density is formally independent of Msusy? since 
5 R cannot interact with sfermions. The Ap bound doesn't allow values for va smaller than 
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FIG. 8: The relic density of a pure higgsino S R as a function of its mass for v§ = 100 GeV, 
va = 1-5 TeV, in the scenario in which the higgsino is the lightest neutralino. The green strip is 
the allowed region by the WMAP at 2a level. 

1.5 TeV, once we assume a relatively small v$- Then the mass of the exchange particle Z R 
becomes heavy, around 1.7 TeV, yielding suppressed cross sections. This is why we obtain 
a large relic density Q x h 2 \ x=s o , which is inversely proportional to thermally averaged cross 
section. Only larger masses for the pure higgsino lead to smaller relic density values. In 
that case, as well as having more phase space available, one is closer to the Z R resonance, 
which increases the annihilation cross section. So, as seen from the figure, the WMAP at 
2<r level is satisfied if the mass of the higgsino lies in 400-450 GeV narrow range. Neither 
b — > s'-f nor (g — 2) M receive contributions from a 5 R higgsino. Unfortunately the WMAP 
allowed region at 2<r level restricts the LSP S R mass to a very small interval. 
We will summarize our result in the next section. 

VI. CONCLUSION 



In this study we have considered the neutralino sector of the LRSUSY model and concen- 
trated on the lightest neutralino state as dark matter, motivated from the fact that WMAP 
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TABLE I: Representative values for the masses of the LSP (Xi)> the next-lightest 
neutralino (\%), lightest chargino (\\ ), and the lightest sfermion (f\) for the three 
scenarios described in this section. The masses are given in GeV. The param- 

eter sets are (M L , M R , My , /x x , /i 2 , M S tjsy , v s ) = (1,1,0.13,1,2,0.2,0.05) TeV for Set 
(1), (M L , M R , My,^, ^ 2 , Msijsy,v S ) = (1,0.2,1,1,4,0.8,0.01) TeV for Set (2), and 
(M l ,M r ,M v ,Hi,H2,Msvsy,v s ) = (2, 2, 2, 2, 0.85,*, 0.1) TeV for Set (3). For all three sets, 
tan/3 = 2, A = 0.2 TeV, and v& = 1.5 TeV are used. "*" indicates that the entry depends 
on MgusY which can be taken arbitrarily large. 



Scenario 


LSP 


M x? 


M x% 


M ± 

Xi 


Ms 

h 


n x h 2 


Parameter Set 


1 


bino 


130.6 


472.2 


500.0 


176.2 


0.106 


Set (1) 


2 


R-wino 


177.1 


924.5 


196.4 


713.4 


0.121 


Set (2) 


3 


5 R higgsino 


412.5 


1426.5 


713.4 


* 


0.115 


Set (3) 



result requires considerations of beyond MSSM models. Now, we would like to summarize 
our findings: 

For the pure states in the LRSUSY model; if the masses of the supersymmetric partners 
(in both bosonic and fermionic sectors) are very small, the lightest neutralino is most likely 
to be a pure B — L bino. This covers regions with M x o ~ 60 — 220 GeV, and Msusy ~ 
150 — 250 GeV. For intermediate to larger values for the masses of the supersymmetric 
partners, the lightest neutralino is most likely to be pure right-wino; this covers regions with 
M x o ~ 150 - 350 GeV, and M SU sy ~ 600 - 1000 GeV. Finally if the mass of the lightest 
neutralino is large, in the 400 — 450 GeV region, the state could be pure (^-higgsino. The 
mass region is severely restricted for pure higgsino masses, but completely independent of 
^susy; as long as Mgusy > 500 GeV, so that none of the sfermions are lighter than the 
lightest neutralino. 

The analysis above shows that the lightest neutralino as pure bino shares some features 
with the similar scenario in MSSM. Both are possible only if the scalars are light, especially 
for My ~ 60 GeV. Outside this region, one obtains a relic density that is too large. This 
similarity is not unexpected: although the bino here is the fermionic partner of the B — L 
gauge boson, and not the hypercharge gauge boson, the same dominant decay channels are 
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open to both binos, so the cross sections are of the same order of magnitude. 

The situation changes when we analyze the case in which the lightest neutralino is the 
right-wino. In MSSM, the case where the lightest neutralino is the left-wino requires a left- 
wino mass of ~ 2.5 TeV to satisfy dark matter limits. For our case, right-wino masses in 
the 150 — 350 GeV regions are in good agreement with the relic density constraints. The 
reason for this manifest difference is that the MSSM left-wino cross section is dominated by 
the decay into W^W^ pairs, which is not available for the right-wino. 

The s-channel is available only for the (5^-higgsino, but this channel is suppressed by the 
large mass of the Zr boson propagator. The cross section is too small for light neutralinos, 
and the relic density is too large. It is less likely that the lightest neutralino would be the 
(^.-higgsino, since only a very small mass range satisfies the dark matter constraints. In 
MSSM, the LSP is even less likely to be the (MSSM) higgsino (see the recent study [32]). In 
LRSUSY, its annihilation cross section into Wl and Z L gauge boson pairs is unsuppressed, 
and the relic density is too small, unless fii > 1 TeV. 

Our scenarios are very different from those present in the NMSSM or in other extensions 
of the MSSM, where a very light bino, or a singlino, or their mixture could be the LSP 
|lllll2l.ll4|. There the cross section is dominated by a light CP-odd Higgs boson or additional 
Higgs resonances, and these particles only couple to no SM particles except for the Higgs 
doublets. 
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APPENDIX: CONSTRAINTS 

We include, for completeness, the constraints on the parameter space coming from the 
anomalous magnetic moment of the muon Aa^ and the branching ratio for b — ► sj. Con- 
straints coming from B s — ► p + p~ are more important than those coming from b — > 57 only 



for tan/5 > 11 |33j, which does not affect the calculation for our choice of parameters. For 
larger tan/? values, we expect the values for the masses of the supersymmetric partners to 
shift to larger values. 

In addition to b — > 57 and Aa M , we have included constraints form the LEP limits on 
supersymmetric masses; we have constrained the mass of the next-to-lightest supersym- 
metric particle to be at least 15 GeV heavier than the lightest neutralino mass (to avoid 
co-annihilations); and we have restricted the Higgs vev's so the p parameter: 



_ x M 2 Zl cos 2 6 W _ 1 cos 2 29 w k\ + k 2 2 



2 1 2 \ 2 

fv-i \~ rCo 

V% + V 2 



(A.1) 



is consistent with its experimental limits p = 0.9998_ 00010 at 2a level j2 



1- (9 ~ 2! 



/' 



In addition to the relic density, another experimental result that can constrain the M1/2 — 
MgusY parameter region is the BNL E821 measurement of the muon anomalous magnetic 
moment a IM = (g — 2)^/2 J34(: 

a M (exp) = 11 659 208(6) x KT 10 (A.2) 

which indicates that the anomalous magnetic moment may need additional contributions 
beyond the SM to be consistent with the experimental values. There are however uncer- 
tainties in how big the deviation from SM really is. A 2.4cr deviation is obtained if the SM 
prediction of hadronic contributions from e + e~ data is used; if instead we use the data from 
indirect haronic r decays the disagreement with the muon anomalous magnetic moment is 
reduced to a 0.9<r deviation. We will use, from the 95% CL. range 35j 

- 0.53 x 10" 10 < Aa M < 44.7 x 10~ 10 (A.3) 
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The dominant contributions to a M in LRSUSY come from the chargino-muon sneutrino and 
neutralino-smuon loops. These are 36J : a^ = a* + a* , where: 

€ = t^'^y.mvM,]^! (A ,4, 

fc=i "*< 

is the chargino contribution. The dominant neutralino contribution coming from neutralino- 
left slepton graphs is: 
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and from the neutralino- right slepton graphs: 



l RH 



7 



rrifj, 



■Myog 



L^i 1(^2 X fc - 
fc=l 

+ 2 



ir, (Ma - 2 tan 2 %iV fe3 ) N k5 ^^ 
V2 m^ 



x 



fe=l ■- ' 

g (JV fcl + tan 2 9 w N k3 ) (N k2 - 2 tan 2 9 w N k3 ) 
> ,(^ - //itan/3) F 4 (y fcAtfl ) l 

— 9 „_9 



9 

Mi 



9 



(A.6) 



M 2 



x fc 

2 ' y k ^L,R 



with the loop functions -F 3 (x), F 4 (z) are given in Eq. (I A. 221) and where x kll 

ml 

M\ 

— 7T- — , and [/, V (and N) are matrices that diagonalize the chargino (and neutralino) mass 
mf, 

Vl,r 

matrices. 



2. b -> S7 



The inclusive decay width for the process 6 — > 57 is given by 



r(6 -> 57) 
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b^F L 
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where M means evolving down to the decay scale \x = m^. The branching ratio can be 
expressed as 

V(h —V c^ 

(A.8) 



£#(& - 57) = V{P ^ Sl) BR SL , 



SL 
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where the semileptonic branching ratio BRsl = BR(b — > ceu) = (10.49 ± 0.46)% and 

r - m b G F\( K CKM)cb\ 2 , x (A Q) 

TsL ~ 192^ M > (A ' 9) 

where z = m 2 c /m1 and g(z) = 1 — 8z + 8z 3 — z 4 — 12z 2 \ogz. The experimental measurement 
from CLEO can be expressed as |38| ]: 

BR(b -> s 7 ) = (3.21 ± 0.43 ± 0.27±°;^) x 10" 4 (A.10) 

The SM contribution is 

AsM = fcn art ( K CKM)ts( K CKM)tb3Xtw [QuFi{x tW ) + F 2 (x tW )} . (A.ll) 

2v 2&fM w 



The matrix elements responsible for the b - 
from the supersymmetric sector of the model 



sj decay acquire the following contributions 



371 ] . For bi decay: 



M lR = Af + A^+A^ a (A. 12) 

with the gluino, chargino and neutralino contributions given by 

Af = -^QMR) E 4-{r&r^F a (^) - ^r^r^F 4 (^ 4 )}, (A.13) 



fc=l d k 
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4?- = -^EEi{( G ^^<)^^*)^M+^^)] 

V2G F i=1 fc=1 m- fc 
+-J« (Gg* - H^){Gf L s - Hf R s )[F,{x Xj ^) + Q u F 4 (* Xj ,J]}, (A.14) 
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yZLrF ■ 1 ; i ''''J 






(V2Gi% R - tf$ L )(V2Gg£ - OF 4 ( V ,)}, (A.15) 



and for the decay of 6#: 

M 7L =^ + ^_+^o (A. 16) 

again, with the following gluino, chargino and neutralino contributions 

4 = -^pGdC(J2) E i{^*^(**4) - ^r» r*F 4 (^)}, (A.17) 

F fc=l dk h 

V2G F j=1 fe=1 m fijk 
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(A.18) 
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where the chargino-quark-squark mixing martices G and if are defined as 

G UL = V jl{^UL)ki, 
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and where the neutralino-quark-squark mixing matrices Gq and H are defined as 
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(A.21; 



y/2m w sm/3 J cos/3 
and where the matrices Tu.d diagonalize the squark mass matrices in the up and down 
sectors. The functions appearing in the expressions above are 

1 



F 1 {x) 
F 2 (x) 
F 3 (x) 
FAx) 



12(x- l) 4 

1 
12(x-l) 4 

1 
2(x-l) 3 



(x — 6x + 3x + 2 + 6x log x) , 
(2x 3 + 3x 2 — 6x + 1 — 6x 2 logx), 
x 2 - Ax + 3 + 2 log x) , 
f x 2 — 1 — 2x log x) 



(A.22) 



2(x-l) 3 

The convention x ab = m 2 a /m 2 is used. C(R) = 4/3 is the quadratic Casimir operator of the 
fundamental representation of SU(3)c- 
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In order to compare the results obtained with experimental branching ratios, QCD correc- 
tions must be taken into account. We assume the renormalization group evolution pattern. 
There is no mixing between left and right-handed contributions. 

A\m h ) = r,- 16/23 {MM w ) + A^(rj^ - 1) + ^ 28/23 - 1)]}, (A.23) 

where r? = a s (m b )/a s (M w ) and A^ = {KcKM)t( K CKM)tb- We choose the 

renormalization scale to be /x = mj = 4.2 GeV. 
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